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We study the two-dimensional toric code Hamiltonian with effective long-range interactions be- 
tween its anyonic excitations induced by coupling the toric code to external fields. It has been shown 
that such interactions allow to increase the lifetime of the stored quantum information arbitrarily by 
making L, the linear size of the memory, larger [Phys. Rev. A 82 022305 (2010)]. We show that for 
these systems the choice of boundary conditions (open boundaries as opposed to periodic boundary 
conditions) is not a mere technicality; the influence of anyons produced at the boundaries becomes 
in fact dominant for large enough L. This influence can be both beneficial or detrimental. In partic- 
ular, we study an effective Hamiltonian proposed in [Phys. Rev. B 83 115415 (2011)] that describes 
repulsion between anyons and anyon holes. For this system, we find a lifetime of the stored quantum 
information that grows exponentially in L 2 for both periodic and open boundary conditions, though 
the exponent in the latter case is found to be less favourable. However, L is upper-bounded through 
the breakdown of the perturbative treatment of the underlying Hamiltonian. 
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I. INTRODUCTION 

An important open problem in quantum information 
concerns the feasibility of a self- correcting quantum mem- 
ory. Finding a system that protects a quantum state 
from decoherence induced by a thermal bath, without the 
need for active monitoring and error-correction, proves 
much more difficult than in the classical case. On the 
most basic level, this is due to the fact that a classical 
bit only needs protection against logical X operations 
while a qubit needs protection against a logical X and 
Z. If a state is stored in a many-qubit system, a desir- 
able feature is topological protection of the stored (qu- 
)bit: we want a logical error X (or Z in the quantum 
case) to necessitate a number of single-qubit errors <j x 
(or cr z ) that scales with L, the linear size of the memory. 
The simplest model that energetically penalizes <j x errors 
and offers topological protection of a stored classical bit 
is the ID ferromagnetic Ising model. The simplest model 
that penalizes cr x and a z errors and topologically protects 
a qubit is given by Kitaev's 2D toric code Hamiltonian 
PQ. In fact, the latter can be mapped exactly to two 
independent copies of the former |2j [3] . Unfortunately, 
both of these systems are not thermally stable. Once a 
topological defect (a pair of domain walls in the ID Ising 
model or a pair of anyons in the 2D toric code) has been 
created, it can spread and lead to a logical error without 
any further energy cost. The lifetime of a qubit stored in 
the degenerate ground states of the 2D toric code is thus 
for any finite temperature upper-bounded by a constant 
independent of L [3 - 5 . The 2D toric code can therefore 
not serve as a 'quantum hard drive'. These difficulties 
can be overcome if the dimensionality of the systems is 
increased. In the 2D Ising model and the 4D toric code, 
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any sequence of single-qubit Pauli operators that leads to 
a logical error has to surpass an energy barrier whose size 
scales with L. Since the number of error paths connect- 
ing two distinct ground states is exponentially large in 
L, these systems are thermally stable below some critical 
temperature T c , meaning that the lifetime of the stored 
information grows exponentially with L [3j EJ [7] . 

Whether a similar degree of protection for a quantum 
state can be achieved in less than four dimensions is not 
clear. One can show that for every 2D local stabilizer 
Hamiltonian the height of the energy barrier separating 
orthogonal states stored in a degenerate ground state is 
upper bounded by a constant independent of L [8HTD], 
ruling out the possibility of using such systems for the 
fault-tolerant storage of quantum information by self- 
correction. In principle, these no-go results leave two 
ways out: Either one abandons the locality of the terms 
in the Hamiltonian or one goes to dimension 3. Indeed, 
both of these routes have been followed in the recent 
literature. While the 3D toric code is not thermally sta- 
ble, Haah showed in a recent breakthrough the existence 
of 3D Hamiltonians with local interactions that have no 
string-like logical operators [11] . Unlike anyons in the 2D 
toric code, defects cannot move further than a certain 
constant distance away without creating other defects. 
This property implies a logarithmically growing energy 
barrier between orthogonal ground states, leading one to 
expect a lifetime that grows polynomially with L [12]. 
However, the best known lower bound on the memory 
lifetime of Haah's Hamiltonian is upper-bounded by a 
constant independent of L [13] and further improvement 
is not expected [14]. Furthermore, a 3D architecture may 
lead to practical difficulties when accessing the physical 
qubits for syndrome measurement and error correction. 

We therefore believe that the most promising route to 
follow in search for a realistic proposal for a quantum 
memory is to start from the 2D toric code Hamiltonian 
and add terms to it that 
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• can be physically motivated, and 

• lead to a memory lifetime that becomes arbitrarily 
large as L — >> oo . 

Long-range repulsive interactions (1/r "-potential with 
< a < 2) between the anyons lead to a logarithmically- 
growing self-consistent mean field gap for anyon creation, 



yielding a polynomially increasing lifetime [5], see Sec. Ill 
below. So rather than trying to find a Hamiltonian with 
a macroscopic energy barrier between orthogonal ground 
states, this approach seeks to suppress the anyon cre- 
ation rate. The toric code Hamiltonian (involving local 
four-qubit couplings) with non-interacting anyons can be 
obtained as an effective Hamiltonian of the Kitaev honey- 
comb model, which involves nearest-neighbor two-qubit 
Ising couplings [15]. How an a = interaction between 
the anyons can be obtained through such a honeycomb 
model coupled to electromagnetic modes has been stud- 
ied in detail in [T6] . see Sec. IV below. 

In an alternative approach it was shown that coupling 
the toric code to a bosonic field leads to an effective grav- 
itational potential between the anyonic defects [TJJ. Be- 
low some critical temperature, all anyons coalesce to a 
single point. However, the time the system needs to ap- 
proach this metastable state and how to best perform er- 
ror correction in this system have not been investigated 
so far. 

Self-correcting quantum memories are usually dis- 
cussed with periodic boundary conditions, giving the 
toric code its name. This way, the complications that 
arise with the possibility of creating unpaired topological 
defects at the boundaries can be avoided. One expects 
that the influence of the boundaries becomes negligible 
if L becomes large enough, which is certainly correct for 
Hamiltonians with local interactions. However, here we 
study memories with long-range interactions between the 
anyons as proposed in [SJ [16] and show that for these sys- 
tems the influence of the boundary becomes in fact dom- 
inant for large enough L. It can be beneficial and detri- 
mental. Specifically, unpaired anyons from the boundary 
lead to an effective bias for anyons from the bulk to move 
to the closest boundary, thus prolonging th e ti me until 
error correction becomes ambiguous (see Sec. III). On the 
other hand, the ability to create unpaired anyons at the 
boundaries halves the energetic gap above the anyonic 
vacuum. This becomes especially relevant if this gap is 
so strong that the anyonic system is basically restricted 
to its ground state and first excited state (see Sec. IV). 

The paper is structured as follows. In Section[TT]we dis- 
cuss how error correction can be performed in the planar 
code (a toric code with open boundaries) in contact with 
a thermal environment. In Section |Hl| we study the influ- 
ence of the boundaries for a Hamiltonian with spatially 
constant repulsion between the anyons, while in Section 
IV an effective Hamiltonian that describes repulsion be- 
tween anyons and anyon holes is investigated. 



II. ERROR CORRECTION IN THE PLANAR 
CODE 

A. The planar code 

A self-correcting quantum memory is supposed to pro- 
tect a quantum state from a thermal environment by 
means of its internal dynamics and without need for ac- 
tive error-correction. A single error correction step may 
be performed before the stored state is read out. We shall 
use here a version of the toric code first introduced in [18] , 
which, contrary to what the name suggests, is not peri- 
odic but does have a boundary. We will refer to this as 
the planar code. Consider a grid with quadratic cells and 
physical qubits placed on the edges, as depicted in Fig- 
ure [l] We call the four qubits around one unit cell a 'pla- 
quette' and the four qubits around a vertex a 'star'. We 
define plaquette operators A p = (cr 2 )® 4 , where the tensor 
product runs over the four qubits around some plaque- 
tte p and star operators B s = (c^)® 4 , where the tensor 
product runs over the four qubits around some vertex s. 
Plaquette operators on the left and right boundary and 
star operators on the top an bottom boundary are tensor 
products of three Pauli operators only. All of these op- 
erators are commuting since they overlap at zero or two 
qubits. Let the space /C C U = (C 2 )® N (N is the total 
number of qubits) be defined as the space stabilized by 
all plaquette and star operators. That is, /Co is the space 
of all states such that for each three- or four-qubit 
plaquette or star operator A p or B s we have A p \i/j) = 
and B s \i/j) = Since all the plaquette and star op- 

erators are independent (unlike in the toric code where 
the product of all plaquette and star stabilizer operators 
is the identity), one easily verifies that dim /Co = 2, in- 
dependent of the height and width of the grid [19 , such 
that one logical qubit can be stored in this space. States 
in /Co are topologically protected. They cannot be distin- 
guished by any local observable and not be evolved into 
each other by any local unitary. 

We then define a Hamiltonian that imposes an en- 
ergy penalty for the violation of every stabilizer condi- 
tion. Let n p = (1 - (a^ 4 )/2, ry = (1 - (a^ 3 )/2, 
n s = (1 - K)^ 4 )/2, and n s> = (1 - {a x )^)/2. The 
tensor products run over the qubits depicted in Figure [I] 
These operators have eigenvalue for states that satisfy 
the corresponding stabilizer conditions and eigenvalue 1 
for states that violate it. For some state |^) G H we 
say that an anyon is present at plaquette p (vertex s) if 
n p \ip) = (n s \ip) = 1^)) an d that no anyon is present 
if n p \ijj) =0 (n s \ijj) = 0), i.e. we interpret stabilizer vio- 
lations as the presence of anyons. We then use the well- 
known toric code Hamiltonian 
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which simply counts the total number of anyons. The 
code subspace /Co, which is the degenerate ground state 
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FIG. 1. A planar code of size L = 8. Depicted are a four-qubit 
plaquette operator n p , a three-qubit plaquette operator n p /, 
a four-qubit star operator n s , and a three-qubit star operator 
n s t. The logical operator X (Z) is given by any chain of 
Pauli operators a x (a z ) connecting the top and bottom (left 
and right) boundary. 



of this Hamiltonian, corresponds to the anyonic vacuum. 
This Hamiltonian is stable against weak local perturba- 
tions in the sense that the lifting of the ground state 
degeneracy through such a perturbation is exponentially 
small in L pQ. 

Starting from the anyonic vacuum, if a qubit suffers 
a spin- flip error <j x (phase- flip error <j z ), two plaquette 
(star) anyons are created on the two adjacent plaquettes 
(vertices). Once an anyon exists, it can move on the sur- 
face without any further energy cost, whereby creating 
further spin- or phase-flip errors along its path. Two 
anyons of the same kind can fuse to the vacuum. On 
the two horizontal (vertical) boundaries of the grid in 
Figure [l] a single plaquette (star) anyon can be created. 
Similarly, the boundaries can absorb single anyons. If 
two anyons are jointly created from the vacuum, move 
around and then fuse to the vacuum, the produced er- 
ror path can be expressed as a product of stabilizers and 
therefore acts trivially on the code subspace /Co- The 
same holds if a single anyon is created on a boundary and 
then is absorbed by the same boundary again. The only 
possibility to act non-trivially on the state stored in /Co 
is if an error path connects the two opposite boundaries. 



We therefore define the logical operators X = (g) a x and 
Z = ® a z acting on all qubits along the paths depicted 
in Figure [l] Since all such products that connect the 
two opposite boundaries are identical up to multiplica- 
tion with stabilizers, the precise form of the path does 
not matter. These operators commute with all the sta- 
bilizers (i.e. they are elements of the centralizer of the 
stabilizer group) though are not products of stabilizers. 
They allow thus to act non-trivially on the state stored 
in /Co without leaving any smoking guns in the form of 
anyons. They satisfy X 2 = Z 2 = I and XZ — —ZX and 
may therefore be seen as Pauli operators acting on the 
encoded qubit. 

At the read-out step, all the operators n s and n p are 
measured and the presence of anyons (the syndrome) is 
detected. The goal is then to annihilate all anyons (ei- 
ther by fusing them to the vacuum or by moving them to 
a boundary that can absorb them) and thereby to undo 
the errors caused by the diffusion of the anyons. More 
precisely, the goal is that the total unitary formed by the 
natural anyon dynamics (creation, diffusion, and annihi- 
lation) plus the error correction procedure is equal to a 
product of stabilizers and thus acts trivially on the code 
space. Even more precisely, the dynamics induced by the 
thermal environment are in fact a probabilistic mixture 
of different unitary evolutions, this will be discussed in 
more detail later on. 



B. Error correction 

We will study planar codes of different sizes L, by 
which we mean that the number of plaquette operators 
from top to bottom as well as the number of star op- 
erators from left to right is given by L. Consequently, 
there are L — 1 four-qubit plaquettes and 2 three-qubit 
plaquettes from left to right and thus a total of L • (L + 1) 
plaquettes. The total number of qubits is then given by 
N = 2L 2 + 2L + 1. It is well-known that in the limit 
of large L error correction in the toric code is unam- 
biguously possible if less than 11% of the physical qubits 
are subject to uncorrelated bit- and phase-flip errors [6]. 
However, qubit errors induced by a thermal environment 
are correlated due to the diffusion of the anyons. Error 
correction therefore typically becomes impossible if a few 
percent of the qubits have suffered errors [5]. 

The error correction step consists of three substeps. In 
a first step, the anyon configuration or syndrome is de- 
termined by measuring all the operators n s and n p . This 
means that we project to the subspaces /Q with a given 
anyon configuration. Since there are L • (L + 1) poten- 
tial locations for anyons of each kind, there are 2 2L '( L+1 ) 
such spaces, each of which has dimension 2 (they can be 
obtained from /Co by applying single-qubit errors), so 
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Indeed, one verifies that 2 • 2 2L "( L+1 ) equals dimH = 
2^ = 2 2L + 2L + 1 . In a second step, a classical computa- 
tion is performed on the error syndrome whose output 
tells how to best annihilate the anyons by fusing them 
with each other or moving them to a boundary of the 
type that can absorb them, which is then done in a third 
step. 

In a more formal language, let po denote the initial 
state stored in the memory with supp po C /Co- The 
influence of errors on this state is then captured by a 
quantum channel (CPTPM) <£ err . The goal of the clas- 
sical computation is then, given knowledge about the er- 
ror model <£>err and the syndrome (the space /Q, that 
is), to find a sequence of single-qubit Pauli operators U{ 
which corrects the errors. In other words, we want U{ 
to map JCi to JCq. Now let Pi denote the projector onto 
/Q. Formally, we can write the error-correction proce- 
dure performed on the corrupted state p t = <£g rr (po) as 
^corr(Pt) = J2i U iPiPtPiUj . At the end of the day, we 
want the error 



5(t) 1= \\P0- ($corrO$* rr ) (p ) || x 



(3) 



to be as small as possible for any given encoded state po. 
In the corrected state (<£>corr o ^ l err ) (po) no anyons are 
left. We therefore either have successfully corrected all 
errors, performed a logical X, a logical Z or both, thus 

($corr O $* rr ) (po) 

= (1 -px)(X -Pz) • Po ~Pz) • Xp X 

+ (1 - Px)pz • ZpoZ + PxPz • XZpoZX (4) 

(assuming that the error correction procedure treats 
plaquette- and star-anyons independently). We therefore 
have 

S(t) < 2 • - Pz) + (1 ~Px)pz + PxPz) 

<2(p x +p z ) . (5) 

In order to obtain simple scalar functions that charac- 
terize the decay of the stored quantum information, we 
study the autocorrelation functions 
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-N 



tr 



X-(^ corr o^ rr ) f (X) 
Z-^ corr o^{Z) 



(6) 
(7) 



The prefactor is such that C^ rr (0) = C c ^ rr (0) = 1, as- 
suming that no operation is performed on the stored in- 
formation if no anyons are measured. 

This assumption is in fact less trivial than it may seem. 
Performing a logical operation in the error correction step 
is beneficial if an odd number of logical operators have 
been performed by the bath. Assume that the bath in- 
duces logical errors that leave no anyons with rate r. 
Then, the probability that no logical error has been per- 
formed is in fact small for times t ^> r _1 . The probability 



that after time t k logical errors have been performed is 
given by the Poisson distribution, 



P(k,rt) = 
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The Poisson distribution is peaked around r£, which may 
be an odd integer. However, the probability 



J2 P(k,rt)= 1 -(l + e-^) 



(9) 



k even 



of an even number of errors is greater than | for any 
rt, such that the optimal strategy is, indeed, not to do 
anything if no anyon is detected. 

C* TT (i) is 1 if after error correction at time t no logi- 
cal Z-operator has been applied and —1 if one has been 
applied (i.e. an odd number of ovoperators has been ap- 
plied to any line connecting the left and right bound- 
ary). Therefore, C* rr (t) = 1 — 2pz and analogously 
C^ orr (t) = 1 — 2px- In conclusion we have 



\Po ~ 

< 



- ($corr 0$*„) (po)^ 

(l-C c ^ rr (t)) + (l-Cf orr (t)) . 



(10) 



We define the lifetime r(e) of the memory as the maximal 
time such that min (C^ rr (t), C^ OTT (t)} > 1 — e for all 
t < t(s), implying that S(t) < 2e for t < r(e). 

The total evolution <£ C orr ° 3>err i s a statistical mixture 
of different unitary evolutions. In the numerical simula- 
tions, we will in each run follow a definite unitary evolu- 
tion, such that C^ orr (t) is at any time given by ±1. The 
probability of a certain unitary is thereby determined by 
the error model ^> e rr- Sampling over a large number of 
runs, we obtain a smooth function C^ orr (i). 

We say that two sequences of single-qubit Pauli oper- 
ators are equivalent if they are identical up to multipli- 
cation with stabilizers. For every given anyon configu- 
ration, there are four equivalence classes of errors that 
produce it from the vacuum. These equivalence classes 
can be mapped onto each other by application of the 
logical operators /, X, Z, XZ. Given a syndrome, the 
goal is to guess the most likely equivalence class of er- 
rors that has produced it, which one allows to remove 
the anyons without disturbing the stored quantum infor- 
mation. Calculating the probabilities of the four equiv- 
alence classes is numerically too costly to be performed 
with current technology. We thus make the simplifying 
assumption that the most likely error path that has led 
to the given syndrome is an element of the most likely 
equivalence class. This may not be true for every possible 
anyon configuration but seems a reasonable approxima- 
tion. Applying stabilizers to the most likely error path 
will produce further error paths with identical or slightly 
lower probabilities that are elements of the same equiv- 
alence class. Rather than finding the error path with 
maximal probability, we may equivalently find the error 
path with minimal weight, if we define the weight to be 
the negative logarithm of the probability that a certain 
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error chain has occurred. Taking the negative logarithm 
ensures that the weight is additive for independent er- 
ror chains. This is known as the Shannon information 
content of an event in classical information theory [20] 
and up to a constant factor the only function having the 
additivity property. 

To illustrate this, let us consider a concrete simple error 
model. We assume that each physical qubit suffers a spin- 
flip error with probability p x and a phase-flip error with 
probability p z , 

cb = $® N 

err single — qubit ' 

^single — qubit 

= (1 -p x )(l -p z ) ■ u+p x (l -p z ) • a x ua x 

+ (1 - Px)Pz ■ cr z ua z + p x p z • (TyU(T y . (11) 

The weight of an error chain involving i x spin-flips and 
i z phase-flips is then 

4 • In + i z . In + const , (12) 

Px Pz 

allowing us to minimize (for p x ,p z < \) t x and t x inde- 
pendently. The number of single-qubit Pauli operators 
necessary to connect two anyons with each other is given 
by the so-called 'Manhattan distance' of the anyons, i.e. 
the sum of the horizontal and the vertical coordinate dif- 
ference of two anyons. Similarly, the weight of a chain 
connecting an anyon to a boundary is given by the hori- 
zontal or vertical distance. In detail, our minimal weight 
matching of n anyons (of one kind) then works as follows. 

1. Perform a Delaunay triangulation on the set of 
anyon coordinates, thereby restricting the full 
graph of n ( n ~ 1 " ) edges between anyons to 0(n) 
edges. Calculate the Manhattan weights of all 
edges in the restricted graph. 

2. For every anyon, add a 'virtual' partner on the 
closer boundary able to absorb it and add an edge 
to the graph with weight given by the distance to 
the boundary. 

3. Connect all virtual anyons with zero- weight edges. 

4. Perform a minimum- weight matching of the graph 
obtained this way. 

Points 1. and 4. are identical to the methods used in 
[5j [21]. We perform the Delaunay triangulation using 
the library Triangle [22] while for the minimal- weight 
perfect matching we employ the library Blossom V [23] 
implementing the 'blossom' algorithm due to Edmond's 
[24] . The numerical cost of this procedure is strongly 
dominated by the last step. Adding the virtual anyons 
ensures that each real anyon can be connected to the 
closest edge able to absorb it and that there is always an 
even number of points in the graph entering the perfect 
matching algorithm. Giving the edges between virtual 
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FIG. 2. The plot shows the autocorrelation function C* OYY 
as a function of the fraction of physical spins p z that have 
suffered spin-flip errors. Circles give the numerical results, 
the lines are guides to the eye. The curves correspond to 
different lattice sizes L = 32 (blue), L = 64 (magenta) and 
L — 128 (red). We anticipate that in the limit L — » oo error 
correction is unambiguously possible for any p z < p c , with 
Pc > 0.102. 

anyons weight zero ensures that those virtual anyons that 
are not connected to a real one can be removed at no cost. 

Employing this algorithm, we obtain the decay of the 
autocorrelation function C^ orr as a function of p x illus- 
trated for different lattice sizes in Figure [2] The curves 
for different lattice sizes intersect for p c > 0.102, so 
that in the limit of L — >> oo error correction is possi- 
ble for p Xl p y < Pc Our numerically obtained value 
is only slightly smaller than the theoretical value p c c± 
0.1094 ± 0.0002 [6]. For the uncorrected or 'bare' auto- 
correlation functions we have 2 _Ar tr [X • <3>| rr (X)] « 
in the whole parameter regime depicted in Figure [2] so 
there is a regime where our error-correction procedure is 
maximally beneficial, bringing the autocorrelation func- 
tion from to 1. 



C. Error model 

We now turn to the situation we are actually physi- 
cally interested in, namely where ^ err is induced by the 
memory Hamiltonian H and coupling of the memory to 
a thermal environment. We consider a Davies weak- 
coupling limit [25] and briefly summarize its discussion 
in [5j[T3j[2T]. In this limit, the evolution of the memory 
is described by a Markovian Master equation 

p t = -i[H, Pt }+jr( Pt ) (13) 

where the interaction between the memory and the bath 
is captured in the unitarity-breaking Lindblad operator 
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C. We assume that the environment is weakly coupled to 
the bath through single-qubit Pauli operators and thus is 
able to induce spin- and phase- flip errors, leading to tran- 
sitions between eigenstates of H that differ only by the 
application of a single-qubit Pauli operator. Processes in 
which an energy uj is transferred from the anyonic sys- 
tem to the bath happen with rate 7(0;), which depends 
on how the bath is modeled. An expression for 7(0;) often 
found in the literature is given by 



7(0;) = 2n n 



1 



-/3a; 



>|/u 



(14) 



and can be derived from a spin-boson model [26j [27] . 
In the following, we set the cutoff frequency of the bath 
uo c — > 00 for simplicity. A bath with n = 1 is called 
'Ohmic', whereas one with n > 2 is called 'super-Ohmic'. 
Only the former case is considered in the numerical sim- 
ulations [28]. To summarize, we will use 



7M = 7(0) 



poo 



1 



(15) 



where we think of 1/7(0) = (2k\T) 1 as the relevant mi- 
croscopic timescale, since the diffusion of anyons is widely 
determined by 7(0). Note that ([14} and ([l5j) fulfill the 
detailed balance condition 7(— Jj = 7(0;) • e _/3u; , guar- 
anteeing that the Gibbs state is the fixed point of the 
Markovian dynamics, C (e - ^) = 0. 

In such a physical model, the weight of a hypothetical 
error chain is not simply given by its Manhattan length, 
as it was the case in (12). It is in general not true that 



the most likely error chain is the one with the smallest 
number of spin flips. The number of hoppings some time 
At after the creation of an anyon or a pair of anyons is 
Poisson-distributed and for large 7(0) • At there is in fact 
a small probability that the number of spin flips is still 
small. Using the Manhattan distance as the weight of 
an error chain connecting two anyons (as done in [5} I2T]) 
seems thus hard to justify. Rather than trying to find the 
most likely error chain, we therefore try to find the most 
likely pairing of the defects (a pairing may either be be- 
tween two anyons or between an anyon and a boundary). 
This is not exactly equivalent to finding the most likely 
equivalence class of errors but should not make a relevant 
difference in practice and is numerically feasible. 

The random walk of an anyon on the grid leads to a 
diffusive spreading of the probability distribution. The 
probability of finding it at time t at a position r relative 
to its position at time t' is 



1 



4irD(t - t') 



2 /AD{t-t') 



(16) 



In the case of an Ohmic bath, the diffusion constant D 
is basically given by the hopping rate 7(0) [5 . Similarly, 
the distance vector of two anyons that have been jointly 
created diffuses with a constant 2D since both of its ends 
are moving. A sensible choice for the weight of an edge 
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FIG. 3. The code has suffered <j x errors at the black qubits in 
a). The syndrome measurement detects plaquette anyons A, 
B, and C. For the error correction procedure, we first perform 
a Delaunay triangulation of the full graph with edges con- 
necting all real anyons (in the case of only three real anyons, 
the full graph and the triangulation coincide). Then, vir- 
tual anyons A 1 , B' , and C' are added. The graph b) is then 
used for the perfect matching algorithm, where the obtained 
matching is highlighted with arrows. Error correction fails, 
because anyons B and C are fused and anyon A is moved to 
the upper boundary, thereby performing a logical X operator 
on the qubit stored in JCq. 



between two anyons is therefore the square of their Eu- 
clidean distance, while for an edge connecting an anyon 
to its closer boundary we take the weight to be twice 
the square of the Euclidean distance. A more thorough 
justification of this choice can be found in Appendix [A] 
An example for the application of the error correction 
algorithm is given in Figure [3] 
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III. LONG-RANGE REPULSION BETWEEN 
ANYONS 

A l/r a repulsive potential with < a < 2 between 
the anyons allows one to increase the lifetime of the toric 
code arbitrarily by increasing L [5, . We study the a = 
case here, since for this case proposals of its physical im- 
plementation exist [5j[T6]. Furthermore, this case can be 
analyzed analytically without need for a mean-field ap- 
proximation and provides the best scaling of the lifetime 
with L. Since the repulsive potential is independent of 
the anyon distances r, its numerical simulation has the 
lowest cost. 

We study a system with a repulsion between the anyons 
which is spatially constant, so the total Hamiltonian is 



and 



H = #i 



Kitaev 



p^p> s ^s' 



where i^Kitaev is as m @ [29 • Basically, the total energy 
is now parabolically rather than linearly increasing in 
the anyon numbers. The physical realization of such an 
interaction and its effect on the lifetime of the toric code 
have been studied in detail in [5j [16j |21] . 

Let us first discuss the toric code, i.e. a memory with 
periodic boundary conditions. The lifetime of the mem- 
ory is given by the time when finding the most likely 
anyon pairing becomes ambiguous. After a time t, the 
distance between the two anyons of a pair is of order 



Dt, with a diffusion constant D. Let n eq be the equi- 
librium density of anyons, such that their average sepa- 
ration is ~ l/y/n eq . We may then estimate the lifetime 
of the memory as the time when the anyons have diffused 
over their average distance and error correction becomes 
ambiguous, 



t(s) 



c(e) 
D • n P , 



(18) 



We use c(e) as a single fitting parameter, which can be 
thought of as a critical fraction of spins affected by errors 
and will be of order of a few percents [5]. 
Let 



A + A(L 2 n eq - 1) 



(19) 



denote energy per anyon in equilibrium. Since there is 
either one or no anyon at each position, the equilibrium 
anyon density can be determined self-consistently from 



[exp(/3e eq ) + 1] 1 



(20) 



The diffusion constant is D = 7(0)+47(— 2e eq ), which for 
an Ohmic bath is widely dominated by the first summand 

Straightforward algebra (c.f. Appendix p| gives two 
simple bounds on n eq for large enough L. We have that 



n eq > J2 if L 2 > exp(/3A) + 1 (21) 



1 

n eq < JJ2Z 



if * A 



L £ - 1 2 - e 
Putting these bounds into ( |20| ) we also find 

d + K L 2 



< exp(pe e 



(22) 



(23) 



for large enough L. The lifetime (|18j) is inverse in n eq and 



will thus for any e > grow faster than L 2 ~ £ as L — >• oo. 

In the case of the planar code, single anyons can be cre- 
ated at and absorbed by the boundaries. The anyon pro- 
duction rate per boundary spin is 7(— e eq ) (if the anyon 
density approaches the equilibrium density), while it is 
27(— 2e eq — A) for spins in the bulk. Again, e eq can be 



determined self-consistently from ([19]) and ([20]), where 
we replace L 2 by L(L + 1) in (19) for the planar code. 



In order to find the total creation rates, these single- spin 
production rates have to be multiplied by 2(L + 1) and 
2L 2 — 1, the number of spins on the boundary and in 
the bulk, respectively. The total rate for production of 
anyons on the boundaries is then, using the bath (|15|), 



2(i + l). 7 (-e eq )~2(L+i; 



/3e e 



exp(/?e eq ) 



7(0) , (24) 



while the total rate for the production of anyons in the 
bulk is 



(2L 2 - 1) • 2 7 (-2e eq - A) ~ (2L 2 - 1] 



4/3e, 



eq 



exp(2/?e eq ) 



7(0) . 
(25) 



Since the anyon density vanishes for large enough L, we 
neglected for these estimates that in fact only those spins 
that have no adjacent anyons should be considered for the 



anyon production rate. Using (23) with e = we see 



that both the anyon production rate on the boundaries 
and in the bulk go to as L — >• oo. But which of the two 



will dominate for large L? Applying again (|23|), it is clear 
that the ratio 



exp(/3e eq ) 



iL of the boundary to the bulk 
rate is in fact diverging with L. The analytics are in ex- 
cellent agreement with numerical simulations (Figure [4| 
that show an increasing ratio after a slight minimum at 
L ~ 16. 

Let us note that the ratio of "active" anyons is in fact 
slightly different than what the above analysis suggests. 
The probability that an anyon created at the bound- 
ary is not reabsorbed before it gets to the second row 



is 



7(0) 



7(0) 



7 (Q)+ 7 (+e eq ) - tT+^) * Similarl y> the Probability that 
a newly created pair is not immediately re-annihilated 
^ s — i(2e^ °+A) ' ^ we are on ^ interested in the ratio 
of boundary to bulk anyons that ever get away from 
their place of creation, we thus obtain a further factor 

7 (2e eq +A) 
67(+e eq ) 

for large enough L. The ratio of "active" anyons created 
at the boundaries to "active" anyons created in the bulk 

is thus finally CTp ff* eq) . 

An anyon created at the boundary has three possible 
fates. 



For the bath (15), this factor converges to 
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FIG. 4. Ratio of the number of anyons created on the bound- 
ary to the number of anyons created in the bulk if the anyon 
population approaches its equilibrium value. The blue circles 
show numerical results sampled over a time 10 5 • OiA)" 1 . 
The red curve shows the analytical prediction exp(/3e eq )/4L. 
We have used p aram eters /3A = 1/0.3 and A/ A = 0.1 and a 
bath 7(cj) as in (15). 



1. It can be reabsorbed by the same boundary (which 
will happen to most of them) or fuse with an anyon 
created at the same boundary, leading to a triv- 
ial operation performed on the degenerate ground 
space. 

2. It can fuse with an anyon created at the opposite 
boundary or be absorbed by the opposite boundary, 
leading to a logical error. 

3. It can fuse with an anyon created in the bulk. 

We will argue that the time it takes to create a logical 
error via possibility 2 is larger than the lifetime of the 
memory in the toric case and that possibility 3 has a 
beneficial effect on the lifetime of the memory. Combin- 
ing these arguments with the fact that anyon production 
on the boundary surpasses the anyon production in the 
bulk, we conclude that the boundary has a positive net 
effect on the lifetime of the memory. 

Let us study the time it takes to create a logical error 
through possibility 2. There are two time-scales involved 
in the creation of such an error: The time it takes to 
create an anyon that walks during its lifetime at least 
once to the opposite half of the grid and the time it takes 
to do this walk. What is the probability that an anyon 
created at a boundary ever moves to the opposite half 
of the grid and what is the expected number of hop- 
pings necessary for this? Up to the boundary conditions, 
this problem is exactly equivalent to a well-known math- 
ematical problem called The Gambler's Ruin. If a newly 
created anyon is to ever reach the opposite half, it first 



has to avoid immediate reabsorbtion. The probability of 
ever getting to the second row is ~ x . Now imag- 

ine that the distance of the anyon to the boundary that 
has created it represents a gambler's bankroll. He starts 
with one unit of money (the anyon starts at the first 
row) and then does a series of fair coin flips for one unit 
of money each. He ends when either he gets broke (the 
anyon gets reabsorbed by the same boundary that has 
created it) or his bankroll reaches L/2 units of money 
(the anyon reaches the opposite half of the grid). Since 
the gambler starts with one unit of money and he does 
a series of games with zero expectation value, his prob- 
ability of ever reaching L/2 has to be 2/L, which also 
answers the corresponding question for the anyon. Note 
that movements of the anyon parallel to the boundary 
that has created it are irrelevant for the production of 
logical errors, which is why the problem can be mapped 
to such a one-dimensional one. Note that in the case of 
anyons this probability provides in fact an upper bound, 
since we neglected the possibility that the anyon fuses 
with another one before reaching the opposite half of the 
grid. A formally correct treatment of the problem can 
be found in Appendix [C] There, we also show that the 
average number of coin flips (hoppings perpendicular to 
the creating boundary) needed to reach L/2 is (L/2) 2 . 
To summarize, anyons are created at the boundaries 



with rate 2{L + 1)7(— e eq ), only a fraction 



7(0) 



do 



eq „ «y - 11 7 (+e eq ) 

not immediately get reabsorbed and only a fraction of at 
most £ ever reaches the opposite half of the grid. If an 
anyon does reach the opposite grid, this takes on average 
a time -2 - (^) (the factor 2 takes hoppings parallel to 
the creating boundary into account). The total lifetime 
of a memory in which anyons are only created at the 
boundaries may therefore be estimated as 



1 



7(+e, 



2(L+l) 7 (-e eq ) 



1 



47(0) 



/3e e , 



7(0) 
f 2L 2 ) 



L 
2 



1 



7(0) 



2 



(26) 



where we used detailed balance. Note that the only 
time-scale that entered the first summand was the anyon 
creation time l/j(—e eq ). However, this cancels exactly 
with 7(+e eq ): the higher the anyon production rate, the 
higher is (by detailed balance) the probability that an 
anyon is immediately reabsorbed, such that the only re- 
maining time-scale is 1/7(0). 

The only property of the bath 7(0;) we used was the 
detailed balance property, so this scaling behavior is in- 
dependent of the particulars of the bath and not spe- 
cific for ( [T5] ) . However, in the case of a super-Ohmic 
bath with 7(0) = the analysis of the Gambler's Ruin 
problem has to be redone with an effective hopping rate 
emerging from indirect hopping processes that scale with 
7(— 2e eq — A) [5 . Since e eq diverges logarithmically with 



L (23), the effective hopping rate becomes vanishing for 



large L, leading to an improved scaling of the memory 
lifetime with L [5]. 
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Using (pCM) we see that the first summand in (26), i.e. 



the timescale needed for the creation of an anyon that 
walks to the opposite half, is for large enough L identical 
to the lifetime of the toric code (18). Furthermore, we 



know from (23) that the second summand grows faster 



than the first one. We conclude that the time needed for 
anyons created on the boundaries to produce a logical er- 
ror grows faster with L than the time before the matching 
of anyons produced in the bulk becomes ambiguous. 

In the limit of very large L, anyons are created almost 
exclusively on the boundary, once the anyon production 
rates have approached their equilibrium values. However, 
as long as the anyon population is small enough anyon 
production in the bulk will outweigh anyon production 
on the boundaries, so that possibility 3 is non-negligible 
even for large L. 

Possibility 3 has the same effect as if the bulk anyon 
with which the boundary anyon fused had moved to and 
been absorbed by the creating boundary. The effect of 
possibility 3 is therefore to cause an effective bias of bulk 
anyons to move to the closer boundary. We saw that 
the probability that a boundary anyon gets a certain 
distance away from its creating boundary decreases at 
least inversely with that distance. This effect is thus the 
stronger, the closer the bulk anyons are to the bound- 
ary. We saw in the discussion of the toric case that error 
correction breaks down if the anyons that have been cre- 
ated as parts of the same pair have moved sufficiently far 
away from each other such that the pair-matching be- 
comes ambiguous. If two such anyons, however distant 
they are, now move to the same boundary this ambiguity 
is resolved without a logical error. 

In conclusion, the equilibrium density of anyons n eq is 
not affected by the boundary conditions, but the fraction 
of anyons created on the boundary becomes dominant for 
large enough L. The time it takes for boundary anyons 
to create a logical error grows faster with L than the 
lifetime of the toric code. We expect a beneficial effect 
from the possibility of boundary anyons fusing with bulk 
anyons. We conclude that for large enough L a planar 
code of size L has a larger lifetime than a toric code of 
size L. 

This is indeed confirmed by our numerical simulations. 
Already for L > 32 the lifetime r(0.1) of the planar code 
exceeds the one of the toric code. For very small mem- 
ories (L < 10) the toric code is superior, which may be 
attributed to logical errors caused by single anyons. In 
the planar code in Figure [I] for example, an anyon cre- 
ated at a boundary needs to perform only 4 hoppings 
to cause a logical error when error correction is per- 
formed. Figure [5] illustrates the temporal decay of the 
stored quantum information by depicting the autocorre- 
lation functions C^ OYY (t) for different lattice sizes and for 
both planar and toric grids. The obtained lifetimes are 
illustrated in Figure [6j 



While we included the analytical prediction (18) for 
the lifetime of the toric code in Figure [6j we cannot give 
such a simple expression for the lifetime of the planar 




FIG. 5. Temporal decay of the autocorrelation functions 
Cf orr (t) for different lattice sizes and boundary conditions, for 
the Hamiltonian (17) and the bath (15). Times are in units 
of (ftiA) -1 and the physical parameters are f3A = 1/0.3 and 
A/A — 0.1. From left to right (at height 0.6, say) we have 
(grid type - size L) toric 32, planar 32, toric 64, planar 64, 
toric 128, planar 128, toric 256, planar 256, toric 512, and pla- 
nar 512. The intersection of C^ rr (t) with the horizontal line 
at height 0.9 is used to determine the lifetimes r(0.1) given 
in Figure [6] 
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FIG. 6. Lifetimes r(0.1) of the quantum information stored 
in the systems described in the caption of Figure [5] The blue 
crosses represent planar codes, the red circles toric codes. The 
sizes L of the codes are 8, 16, 32, 64, 128, 256, 512, and 
1024. The red line shows the analytical prediction for the 
lifetime in the toric case obtained from ( |18| ) with c(0.1) = 
5.1%. [5, find with an analogous plot [Fig~6] c(0.1) = 4.4%. 
The increased lifetime is due the choice of the square of the 
Euclidean distance rather than the Manhattan length as the 
weight of an error chain. 
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code. When L is increased, the effect of the boundary 
changes from adversarial to beneficial. 

IV. THE HONEYCOMB MODEL AS A 
QUANTUM MEMORY 

A. Prom the honeycomb to the planar code 

The toric code Hamiltonian (involving local four- 
qubit couplings) can be realized as a low-energy effective 
Hamiltonian of the Kitaev honeycomb model, which in- 
volves only local two-qubit couplings [15j [16l [31] . The 
Hamiltonian of the honeycomb model can be simplified 
through a spin to hard-core boson transformation [31]. 
In this new language, the toric code Hamiltonian (HI) 
emerges as a low-energy (no hard-core bosons present) 
fourth-order effective Hamiltonian. The qubits of the 
toric code obtained this way are not identical to the phys- 
ical qubits of the underlying honeycomb lattice, but are 
effective qubits. The interactions of ([!]) between four of 
these effective qubits are mediated through processes in 
which two pairs of virtual hard-core bosons are created 
and then fuse again to the vacuum. An interaction in- 
volving three effective qubits would correspond to a pro- 
cess in which a pair of virtual hard-core bosons is created 
from the vacuum, followed by a hopping and an annihila- 
tion process. Clearly, such processes give no contribution 
to the effective Hamiltonian and the third-order effective 
Hamiltonian vanishes. There is no way of obtaining from 
the honeycomb model the three-qubit boundary terms 
introduced in [18] and used so far in this paper. 

Figure [7] schematically shows a planar code obtained 
from the honeycomb model, involving only four- but no 
three-qubit operators. The ground space of the Hamil- 
tonian obtained this way, i.e. the space stabilized by all 
four-qubit operators (c^) 04 , (c^)® 4 , has then a large de- 
generacy growing with the size of the memory, even if we 
forget about the four isolated qubits in the corners that 
do not interact with anything. 

Let Q denote the group generated by all four-qubit sta- 
bilizers. In an abstract language, elements of the Pauli 
group on all effective qubits (except the four isolated ones 
in the corners) that are in the centralizer C(Q) of the sta- 
bilizer group Q, but not in Q itself, perform logical opera- 
tions on the stabilized subspace that cannot be detected 
through violation of a stabilizer. The goal is then to find 
Pauli-like observables X, Z G C(Q) \ Q which allow us to 
store a qubit that is topologically protected. Elements of 
C{Q) \ Q that commute with X and Z will not do any 
harm to the stored qubit, but those that do not commute 
with X or Z can. For topological protection we therefore 
want 

• X and Z to have a macroscopic distance (i.e. to 
necessitate O(L) single-qubit operations), and 

• all elements of C{Q) \ C((X,Z}) to have macro- 
scopic distance. 



2 




FIG. 7. An effective planar code obtained from a honeycomb 
lattice with a boundary. The black dots represent effective 
qubits emerging from the underlying honeycomb model. The 
four isolated effective qubits in the corners are depicted for 
completeness but not needed for information storage. Blue 
(dark) squares (like p) are plaquette operators, red (light) 
squares (like s) are star operators. There are two non- 
equivalent ways of defining a pair of macroscopic Pauli-like ob- 
servables, the "big" logical operators X and Z and the "small" 
logical operators X' and Z' . However, an undetectable log- 
ical error consisting of only three single-qubit errors can be 
performed on all of these operators. For instance, a sequence 
of a x errors on the three qubits along path 1 (2) causes an 
error on the logical Z' (Z) operator. 



Unfortunately, in the planar code obtained from the hon- 
eycomb model there is no choice of logical operators X 
and Z that fulfill these requirements. As illustrated in 
Figure the pairs X, Z and X',Z' fulfill the first re- 
quirement but not the second one. 

We will therefore expand Q in such a way that all ele- 
ments of C(Q) \Q have a macroscopic distance. For the 
space stabilized by Q being non-trivial we need Q to be 
Abelian and to not contain —I. All of these requirements 
are met if we add to Q the single-qubit operators a x act- 
ing on all qubits surrounded by a red circle in Figure [7] 
and the operators <j z acting on all qubits surrounded by 
a blue circle. Error paths 1 and 2 in Figure [7] now no 
longer are elements of C(Q) \Q: 1 is no longer an element 
of C(Q) since it anti-commutes with two blue single-qubit 
operators in Q while 2 becomes an element of Q (it is the 
product of the two red operators enclosed by it). There- 
fore, the "small" operators X' and Z' now satisfy our 
requirements for topological protection while the "big" 
operators X and Z still violate the second requirement. 

Our new stabilized subspace topologically protects ex- 
actly one qubit. In fact, our new code is exactly equiv- 
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alent to the planar code with three-qubit operators on 
the boundary, since the single-qubit stabilizer operators 
effectively eliminate the two degrees of freedom of the 
qubits on the boundary. However, the Hamiltonian dy- 
namics are different since there is no energy penalty as- 
sociated with the violation of the single-qubit stabiliz- 
ers on the boundary. We therefore assume that we are 
able to perform the measurements corresponding to the 
single-qubit stabilizers at the read-out step. Note that 
these operators act on effective qubits, but a x and a z 
measurements performed on them indeed correspond to 
measurements of a single physical qubit of the underlying 
honeycomb lattice. We refer to [16] for details about the 
mapping between physical and effective qubits. 

If one of the single-qubit stabilizers on the boundary 
is violated, we may imagine that an anyon is present at 
the position adjacent to it which is outside of the actual 
grid. For example, a a x (a z ) error can create two pla- 
quette (star) anyons at positions p and p' (s and s') in 
Figure U\ An anyon at position p' (V) can only escape to 
p (s). In terms of anyon dynamics there are now for both 
kinds of anyons two boundaries that can create and ab- 
sorb them and two boundaries that "attract" and store 
them. For example, a plaquette (star) anyon at position 
p (s) can reduce its energy by hopping to position p' (s f ). 
Analogously, escaping from one of these boundary posi- 
tions has an energy cost. An anyon that hops to position 
p" or s" has been absorbed by the corresponding bound- 
ary and can, unlike one at position p' or s' no longer be 
detected through a violated stabilizer operator. During 
error correction, an anyon stored at position p' or s' will 
be moved to the interior of the grid. 



B. Repulsion between anyons and anyon holes 

As discussed, the simple toric code Hamiltonian 
does not provide a lifetime of the stored quantum in- 
formation that can be increased by making the memory 
larger. Pedrocchi et al. [16] studied a honeycomb model 
as introduced above coupled to cavity modes. The cavity 
modes allow the read-out of the error syndrome of the ef- 
fective toric code through frequency shifts. In a resonant 
parameter regime, the Hamiltonian 

# eff = A^n a rv (27) 

a, a/ 

is found perturbatively [eq. (52) in [IB]]. Here, the sum 
runs over all stars and plaquettes (four-qubit oper- 
ators only) and n a = 1 — n a counts anyon holes. For 
A > this Hamiltonian describes an effective repulsion 
between anyons and anyon holes. For details about how 
A > can be achieved, we refer to [IB] and Section VII 
in 0. 

The requirement to stay strictly in the perturbative 
regime puts an upper bound L* on the linear size L of 
the memory [IB] . 



With Hamiltonian ( |27| ) and a total of N four-qubit op- 
erators, the gap above the anyonic vacuum is (N — 1)A 
and as long as n a < y the cost to add a further anyon 
is at least |A. If L denotes the number of stars from 
top to bottom and the number of plaquettes from left to 
right (so L = 3 in Figure [7) we have N = 2L(L + 1). 
With an anyon creation gap that grows quadratically in 
L and a natural bath like ( [15] ) in which the creation rate 
decreases exponentially with the gap, we expect the pro- 
duction of anyon pairs in the bulk to be negligible against 
the creation of anyons at the boundaries. Note that if 
we say that a plaquette anyon is created at a bound- 
ary where single-qubit a z measurements are performed, 
this implies that furthermore a virtual anyon is "stored" 
in this boundary. From the discussion of the Gambler's 
Ruin problem in Appendix[C]it is clear that an anyon will 
do at most ~ L 2 hoppings before being absorbed by or 
stored in a boundary. We thus expect the time an anyon 
is present (and not stored in a boundary) to be much 
smaller than the time it takes to create a single anyon. 
To summarize, we are already for moderate values of L 
in the regime 

^ « [4L • 7 ((N - lJA)]" 1 « [L 2 • 7 (2(JV - 2)A)] " X 

(28) 

where almost always one or no anyon is present. 
In this regime, error correction works as follows. 

1. Move any anyons in the interior of the grid per- 
pendicularly to the closest of the four boundaries. 
(Each anyon is thereby moved to its most likely 
place of creation. We never assume that a pair of 
anyons has been created in the bulk.) 

2. If an anyon is "stored" in a boundary (at position 
p' or s' in Figure [7[ say), move it to the interior of 
the grid (to position p or s). 

3. The anyons at each boundary can be matched with 
each other and the two- adjacent boundaries in two 
different not trivially suboptimal ways. Chose the 
one with lower weight, where the weight is given by 
the square of the Euclidean distance. 

For this system, the perfect matching problem could be 
reduced to a one-dimensional one. Both possible match- 
ings that are not trivially suboptimal can be explic- 
itly tested, such that no approximative algorithms are 
needed. 

Adding single-qubit measurements at the boundary 
spins ensures that in a planar code without three-body in- 
teractions all error paths that lead to a logical error after 
error correction consist of O(L) single-qubit errors. Fig- 
ure [8] illustrates schematically all error paths of a single 
anyon that do lead to a logical error after error correction 
and all that do not. Strictly speaking, this figure is only 
valid for the first anyon, as a concatenation of paths that 
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FIG. 8. A planar code consisting only of four- and one-qubit 
stabilizers. We only consider a x -errors and the corresponding 
plaquette operators here. As in Figure [7] a z -measurements 
are performed on the top and bottom boundary. The top and 
bottom boundary are thus able to store and detect the pres- 
ence of plaquette anyons. The logical Z operator is a path 
of single-qubit a z operators connecting the top and bottom 
boundary (c.f. the operator Z' in Figure [7]). The dashed line 
lies midway between the left and right boundary. All error 
paths in the left part of the figure will not lead to a logical er- 
ror after error correction, while for the error paths in the right 
part of the figure a logical X operator is performed. Whether 
an error path of type (*) is correctable depends not only on 
its length, but also on its position relative to the boundaries. 
However, any error path of type (*) and of horizontal length 
below L/2 is correctable. 



do not cause an error may lead to an error after error 
correction. 

It is in principal possible to find an analytical expres- 
sion for the probability that the first anyon causes an 
error in such a planar code of size L. However, we ex- 
pect from the discussion of the Gambler's Ruin prob- 
lem in Appendix [C] that this probability takes the form 
7 (o) +7 7 ((Ar_i)A) ' Tji u P t° some constant factor of order 
O(l). This probability will, due to both factors, be small 
such that, following the discussion before (10), we expect 




-0.2 
10° 



the error-corrected autocorrelation function after the cre- 
ation of n single anyons to take the form 



C^ OYY (n) 



7(0) 



L 7 (0) + 7 {{N -1) A) 



(29) 



where c' is a constant. Of course, this approximation can 
only be valid as long as C^ orr (n) is still relatively close to 
1. The numerics in Figure [9] are in excellent agreement 
with this prediction as long as C^ orr (n) > 0.5. 

Correspondingly, the number of single anyons that are 
created during the lifetime r(e) of the stored quantum 
information will be 



e-L 7 (0) + 7 ((TV - 1) A) 
d ' 7 (0) 



(30) 



This prediction is compared with numerical simulations 
in Figure [To] Since there are for each kind of anyon 
2 • (2L + 1) qubits where a single anyon can be created 



10 1 io 2 io 3 

number n of single anyons 



FIG. 9. A planar code of size L — 64 consisting only of 
four- and one-qubit stabilizers in contact with a bath with 
7 ((N - 1)A) /7(0) = 40, where N = 2L(L + 1). The verti- 
cal axis shows the number n of single anyons that have been 
created. The green curve shows the error-corrected autocor- 
relation function C^ rr (n) and the blue curve the uncorrected 
autocorrelation function (Q without <3> C orr). Both curves are 
sampled over 1.2 • 10 4 experiments. The red curve shows the 
analytical prediction ( |29[ ) with c = 1.12. We see that per- 
forming a single error correction step before the read-out of 
the stored quantum information allows to enhance its lifetime 
(r(0.1), say) by more than an order of magnitude. 



10 4 
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size (L) 



FIG. 10. The number of single anyons created during the 
lifetime r(e) of the memory as a function of the lattice size L. 
The points show numerical results obtained for a bath with 
7 ((N - 1)A) /7(0) = 40, where N = 2L(L+1), the lines show 



the analytical prediction (30) with c' — 1.12. From bottom 
to top we have e = 0.1, 0.2, 0.3. 
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(including the ones on which single-qubit measurements 
are performed), we obtain a lifetime 



r(e) 



e-L 7(0)+7((JV-l)A) 

d ' 7(0) 
•[2.(2i + l). 7 (-(iV-l)A)]- 1 

£ exp ((N - 1)/3A) 
: 4d 7(0) ' 



(31) 



In the last step, we applied detailed balance and 7(0) <C 
7 ((TV - 1)A). Since N = 2L(L + 1), the lifetime of the 
memory grows exponentially in L 2 . This provides, to 
the best of our knowledge, the best scaling of a quantum 
memory lifetime with the linear size found so far in at 
most three dimensions. However, we have noted that 
strictly speaking L is bounded through the breakdown of 
the perturbative treatment of the underlying honeycomb 
Hamiltonian. So similarly as in [13] the found scaling of 
the lifetime is only valid up to some optimal L*. Still, 
our optimal L* is generic in the sense that changing it 
by a small integer will not have a drastic effect on the 
lifetime. 



We note that the only time-scale that entered (31 ) was 
the time needed to create a single anyon, while 7(0) only 
entered through a probability. The former canceled by 
detailed balance such that the only remaining time-scale 
was the hopping time 1/7(0). In the toric case, anyons 
can only be created pairwise, such that the gap above the 
anyonic vacuum is increased from (N — 1)A (as in the 
planar case) to 2(N — 2)A, leading to a lifetime increas- 
ing as ~ exp(2(7V — 2)/? A), in contrast to ( [31] ). Besides 
complicating error correction, the realistic case of open 
boundaries thus also reduces the exponentially increas- 
ing factor in the lifetime of the memory to (almost) its 
square root. 



C. Topological order at finite temperature 

The above results apply to the code during its ther- 
malization. However, it is also interesting to study its 
properties once it reaches thermal equilibrium. Clearly 
the memory will have completely decohered by this point, 
since the probability of the system being in any of its 
ground states will be equal. However, we can assess 
whether topological order is present. It is known that, for 
the non-interacting case of ([!]), the thermal state is not 
topologically ordered for any finite temperature. How- 
ever, here we show that the interactions of (27) allow the 



topological order to remain stable for all finite tempera- 
tures at which the perturbative derivation of the Hamil- 
tonian (27) is valid. 



To determine whether topological order is present we 
can use one of the topological order parameters designed 
for mixed states, such as the anyonic topological entropy 
[32] . This requires the plaquettes and vertices of the 
code to be split into three regions, A, £?, and C. These 
can be defined arbitrarily, except that A and B must be 



bounded, B must enclose A, C must enclose B and the 
number of plaquettes and vertices in each region must be 
0(L 2 ). The value of the entropy depends on how well 
the anyon configuration within B can be used to deduce 
the net anyonic occupation of A. If this can be done 
perfectly, the entropy takes its maximum possible value 
(for the planar code) of 2 In 2, signaling that the state is 
topologically ordered. 

Let us consider the thermal state of the Hamiltonian 



(27). The Hamiltonian is symmetric under exchange of 



anyons and anyon holes. Since the gap above the ground 
states with all anyons or all holes grows as L 2 , the ther- 
mal state in the thermodynamic limit will be an equally 
weighted mixture of the all hole or all anyon states. As 
such, measurement of the occupancy of any plaquette can 
be used to determine the occupancies of all plaquettes. 
The region B used in the definition of the anyonic topo- 
logical entropy can then clearly deduce the net occupancy 
of the region A. The value for this entropy therefore takes 
its maximum value of 2 In 2, and the state is found to be 
topologically ordered. 



V. CONCLUSIONS 

Stabilizer Hamiltonians with local interactions in 2D 
do not, and in 3D seem not to, offer the possibility 
to passively store quantum states for a time that can 
be made arbitrarily larger than the relevant microscopic 
time-scales. Inducing long-range interactions between 
the excitations of a 2D Hamiltonian whose ground states 
are topologically ordered seems thus the most promis- 
ing approach towards a realistic proposal for a quantum 
memory. For such long-range interactions, the influence 
of the boundary of the memory (which every realistic 
memory will have) is not negligible even for large L. We 
discussed two Hamiltonians proposed in the recent liter- 
ature and showed that for those the boundary becomes 
in fact dominant. Operationally, this fact becomes rele- 
vant during the error correction step before the read-out 
of the stored quantum state. We showed that the clas- 
sical algorithm that determines how to best remove the 
anyonic defects has to depend on the error model, the 
memory Hamiltonian and the boundary conditions. 

With long-range repulsion between the anyons, the en- 
ergy to add a further anyon increases with the number 
of already existing anyons, leading to a vanishing anyon 
production rate. If the production rates approach their 
equilibrium values, the production of unpaired anyons 
on the boundaries will outweigh the production of anyon 
pairs in the bulk. The influence of the boundary anyons 
is beneficial since they lead to an effective bias of the bulk 
anyons to move towards the closer boundary. 

We discussed how a planar code with topological pro- 
tection of the stored qubit can be obtained from a honey- 
comb model with two-qubit Ising coupling and the abil- 
ity to perform single-qubit measurements on boundary 
qubits. In a resonant regime of a coupling of the honey- 
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comb to cavity modes a very strong suppression of the 
anyon creation rate is obtained. Furthermore, most cre- 
ated anyons will immediately be reabsorbed by their cre- 
ating boundary. In conclusion, we have found in this 
regime a lifetime that grows exponentially in L 2 , allow- 
ing in principle to reach macroscopic storage times even 
at moderately high temperatures. The non-local anyon 
interactions in the obtained effective Hamiltonian are so 
strong, that the system is topologically ordered at any 



finite temperature for large enough L, as long as the per- 
turbatively derived Hamiltonian (27) describes the dy- 
namics of the memory accurately. 
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Appendix A: Determining the weight of the edges 
when anyons perform a random walk 

Given an anyon at time t at position a — {a^a^)^ 
where a± is the distance from the left boundary and a2 the 
distance from the upper boundary, what is the probabil- 
ity Pr [a, t' ', t] that it has been created at a time < t f < t 
at the upper boundary? The answer can be obtained by 
summing ( [16] ) over all spins 0, . . . , L on the upper bound- 
ary, so 



Pr [a, t',t] 

rL+l 



I 

JO 



dx 



1 



4ttL>(£ - t') 



-(al + (a 1 -x) 2 )/4D(t-t') 



4 y/^D(t - V] 



■ e 



Erf 



2y/D(t-f) 



-a%/iD(t-t') 



Erf 



L+l-ai 

2y/D(t-t') 



(Al) 



The total probability that that an anyon that is at time 
t at position a has been created on the upper boundary 
is then 



Pr[a,t] = / dt' Pr[a,t',t] • 7 bo 
Jo 



ridary 



(A2) 



where 7boundary(£ / ) is the creation rate of anyons on the 
boundary at time t' and £(t',t) is the probability that 
an anyon that has been created at a boundary at time 
t' does still exist at time t. Unfortunately, the time in- 
tegration cannot be performed in closed form even if we 
take 7boundary and £ to be constant and {...} = 2 in (Al ) 
(corresponding to the L — >• oo limit). 

Still, the above analytics is enough to find reason- 
able weights for the edges in our error correction algo- 
rithm. First, let us note that we can find a similar ex- 
pression like the one above for the situation where we 
have anyons at postion a and b and are interested in 
the probability that they have been jointly created in 
the bulk. The exponential factor becomes in this case 

~ exp — (a — b) 2 /SD(t — t') . Intuitively, the distance 

vector of the two anyons performs a diffusive motion with 
diffusion constant 2D. A sensible choice for the weight 
of an edge between two anyons is therefore the square 
of their Euclidean distance, while for an edge connecting 
an anyon to its closer boundary we take the weight to be 
twice the square of the distance. All terms depending on 
time, the position of the anyons relative to the bound- 
aries and the creation rates give then additive logarithmic 
correction terms to these weights. 

Note that in the matching which is obtained at the end 
of the algortihm the number of edges connecting two real 
anyons is always equal to the number of edges connecting 
two virtual ones. Furthermore, adding the same term to 
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the weights of all edges or multiplying all weights with 
the same positive term does not change the result of the 
algorithm. We may therefore give the edges that connect 
two virtual anyons a weight which takes all logarithmic 
correction terms into account. 

Formally, let d 2 denote the square of the Euclidean dis- 
tance of the pair we are interested in. Then, giving edges 
that connect two real anyons a weight a • d 2 + /?2 (with 
a > 0), edges that connect a real anyon with a virtual 
one a weight a • 2d 2 + j3\ and edges that connect two 
virtual anyons a weight is equivalent to giving them 
weights d 2 + (/?2 — 2d 2 and —j3\/a respectively, 

which is again equivalent to giving them weights <i 2 , 2d 2 
and (/?2 — 2/?i)/a respectively. Rather than calculating 
the term (/?2 — 2/?i)/a analytically (which we cannot) 
we may then consider the weight of the virtual edges as 
a single optimization parameter of our algorithm. We 
have found numerically that in the systems we are inter- 
ested in varying the weight of the purely virtual edges 
offers hardly room for improvement of the memory life- 
time. As before, we will thus take their weight to be zero 



throughout Sec. Ill However, there are certainly regimes 
where these weights become relevant. If, for example, the 
rate for creation of anyon pairs in the bulk is vanishing 
against the rate for creation on the boundaries, the edges 
between virtual anyons should be given a large weight. 



Appendix B: Bounds on the self-consistent anyon 
density 



Using (19) to eliminate e eq in (|20j) we find 



1 = n eq • [exp (f3(A + A(L 2 n eq - 1))) + l] := /(n eq ) . 

(Bl) 

Since f(x) is monotonically increasing in x we have that 
x < n eq 4=> f(x) < 1. We have 



/(-l) = -l[exp(/3A) + l] 



(B2) 



and 



f(j^) > ■ exp (f3A(L £ ~ 1)) , (B3) 



yielding the estimates (21) and (22). 



is biased to move towards the boundary if it is at a posi- 
tion next to it. This is in particular the case immediately 
after the anyon has been created. 

We therefore study the following problem. Consider a 
lattice with rows 0, . . . , L — 1 and an anyon initially in 
row (it has just been created). We are only interested 
in hoppings between the rows and ignore all movements 
within the same row. In row there is a probability p = 
7 (absorb) +7(0) the anyon is absorbed by the adjacent 
boundary and a probability 1 — p that it hops to row 1 
(similarly in row L — 1, but this will be irrelevant here). 
In all other rows, the anyon hops to both adjacent rows 
with equal probability. What is the probability that the 
anyon reaches row L/2 (let L be even for simplicity) at 
least once before it is reabsorbed? 

Let Pi the denote the probability to reach row L/2 from 
row i. The boundary conditions are 



Po = (l-p)-pi 



and 



PL/2 



1 . 



(CI) 



(C2) 



For 1 < i < L/2 — 1 we have Pi — \(pi_\ + or 
equivalently pi+i — Pi — Pi — Pi-i- We find 



L/2-1 

1 = Pl/2 = 



Pi) +Pi 



(L/2-l)>(p 1 -p )+p 1 



(C3) 



Using (CI) to eliminate p\ we arrive at 

» - (Z+TT+T ' (C4) 

which for 1 — p <C 1 (i.e. 7(0) <C 7(absorb)) simplifies to 

2 4 7 (0) 2 



Po ^ (1 -p) • 



L 7(absorb) L 



(C5) 



Now let rti denote the expected number of hoppings 
(perpendicular to the boundary under interest) necessary 
to reach row L/2, assuming that it is eventually reached. 
The boundary conditions are 



n 



(C6) 



and 



Appendix C: The Gambler's Ruin 

In the usual version of the Gambler's Ruin problem, 
the probability of winning or losing is the same for every 
amount of money the gambler possesses. In the anyonic 
hopping problem we are interested in, the probability of 
moving towards the boundary is not excactly indepen- 
dent of the distance to the boundary. Since the rate for 
absorbtion of an anyon by a boundary 7 (absorb) is usu- 
ally much higher than the hopping rate 7(0), the anyon 



n L/2 = . (C7) 

For 1 < i < L/2 — lwe have the recursion = 1 + 
^(rii-i +77,3+1) or equivalently n^+i — rti = ni — ni-i — 2. 
From this we expect a quadratic expression for Ui and 
thus make the Ansatz 

ni=a-i 2 + b-i + c (C8) 

with the unique solution a = — 1, b = and c = (L/2) 2 . 
We conclude that no = (L/2) 2 . 



